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We investigate theoretically the four-wave mixing of optical and matter waves resulting from the 
scattering of a short light pulse off an atomic Bose-Einstein condensate, as recently demonstrated 
by D. Schneble et al. [ Science 300, 475 (2003)]. We show that atomic "pair production" from 
the condensate results in the generation of both forward- and backward-propagating matter waves. 
These waves are characterized by different phase-matching conditions, resulting in different angular 
distributions and temporal evolutions. 
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The interaction between optical fields and atomic 
Bose-Einstein condensates has attracted much recent at- 
tention, due of its importance in the preparation, manip- 
ulation, and detection of condensates, as well as because 
of its interest in fundamental studies of the nonlinear in- 
teraction between Maxwell and Schrodinger waves. A 
number of phenomena have already been studied both 
theoretically and experimentally, including matter-wave 
superradiance 0, Q> tl > coherent matter- wave amplifica- 
tion U IE E| and Bragg spectroscopy 0, @, @ • 

In a trailblazing experiment 0] , a cw laser beam shined 
on a cigar-shaped condensate resulted in the generation 
of a fan-like pattern of momentum sidemodes of the con- 
densate. The initiation of this pattern can be understood 
in terms of a four-wave mixing process involving two op- 
tical fields - the laser field and a so-called end-fire mode; 
and two matter- wave modes - the condensate and a mode 
of momentum such that energy-momentum conservation 
(or phase matching) is satisfied. We recall here that the 
end-fire mode, first predicted by Dicke 1.0] in his study of 
superradiance, corresponds to the privileged direction for 
spontaneous emission along the long axis of the conden- 
sate. Because of momentum conservation, it is clear that 
the generated sidemodes must be in the forward direc- 
tion at a 45° angle between the direction of the incident 
laser and the long axis of the condensate. The subsequent 
generation of further sidemodes results simply from wave 
mixing involving an already excited sidemodc instead of 
the initial condensate at rest. 

Recently, that same MIT group ^1 reported the re- 
sults of an experiment where the incident cw laser was 
replaced by an optical pulse. This led to the remarkable 
result that for short enough pulses, backward-scattered 
atoms (with a momentum component antiparallel to the 
direction of the pump field) were also observed. More- 
over, the backward peaks exhibited a slightly different an- 
gular distribution compared with the forward peaks. The 
qualitative difference in diffraction patterns for the short- 
and long-pulse regimes was attributed to the transition 
from the Raman-Nath to the Bragg regime of diffraction, 
that is, to the onset of energy-momentum conservation 
for long enough interaction times. 

In this Letter, we present a theoretical description of 
this experiment based on an extension of the quasi-mode 



approach of Ref . 0] . We give a full dynamical treatment 
of both the vacuum photon modes and the condensate 
sidemodes that interprets the diffraction pattern in terms 
of atom-photon wave mixing and show explicitly that 
the counter-propagating (backward and forward) matter- 
wave quasi-modes result from the quantum-correlated 
parametric excitation of atomic pairs. The difference 
in their angular distribution and dynamics results from 
their distinct phase-matching conditions. 

A general theoretical framework to describe the inter- 
action of ultracold atoms with light waves was presented 
by Zhang and Walls in Ref. ■ I* 1 the situation at hand, 
the atomic system, assumed to be at zero temperature, 
interacts with a far off-resonant classical laser field of 
(real) Rabi frequency fiz,, wave vector and frequency 
u)l = cltL, as well as with a continuum of electromagnetic 
field modes of wave vector k and polarization A charac- 
terized by the bosonic annihilation operators -BkA- In the 
case of large atom-laser detunings A, we are justified in 
adiabatically eliminating the excited atomic levels, leav- 
ing us with just the bosonic ground-state matter-wave 
field operator -0i(r,t). We furthermore perform the ro- 
tating wave approximation and express the dynamics of 
the coupled atoms-radiation system in a frame rotating 
at the pump laser frequency to find 
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where = i\/^^>k / {KV)d-e^\ is the coupling strength 
of the atom with the corresponding vacuum mode (V 
is the quantization volume, e^A the photon polarization 
unit vector and d the atomic dipole moment). In Eq. (1), 
the atomic Hamiltonian Ho(r) includes the usual kinetic 
and trapping potential terms, the ac Stark shift arising 
from the pump laser as well as nonlinear atom-atom in- 
teractions resulting e.g. from two-body collisions. 

In the absence of electromagnetic fields, the conden- 
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sate at zero temperature is taken to be in the ground 
state Vo( r ) 01 -^07 with (H — %^,)<p (r) = 0, being the 
chemical potential. The interaction of this condensate 
with light shifts the atomic momentum by the photon 
recoil, and transfers condensate atoms to states of the 
form (fio(r) exp(iq-r). For condensates large compared to 
an optical wavelength, these states are still approximate 
eigenstates of Ho, with eigenenergy shifted by the recoil 
frequency cu q — hq 2 /2M, i.e., H [(fio{ r ) exp(iq • r)] ss 
h(fi + Lo q ) Vo( r ) exp(jq • r), an approximation similar to 



the slowly varying envelope approximation in optics. Fol- 
lowing Ref. 0| , this suggests expanding tpi (r, t) as 

Mr, t) = e-^ PoWe^'ci,, (3) 
q 

where the atomic quasi-mode field operators approxi- 
mately obey bosonic commutation relations, [c q ,cL] = 
<5q,q'- Inserting this expansion into Eq. Q gives 
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where n q (k) = J d 3 r |<^(r)| 2 exp[— i(q+ k — k^) • r] is the spatial Fourier transform of the condensate density profile. 
Consistently with the slowly- varying approximation implicit in the expansion JSJl , it satisfies the approximate relation 
n*(k)n q '(k) w |n q (k)| 2 <5 qiq <, a condition that indicates that the discretization of q is such that there is neglibible 
overlap between shifted momentum wave functions. This property is used repeatedly in the following. 
We now proceed by formally integrating Eq. (J2J to find 
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where we have used Eq. J2J and have adopted the Markov 
approximation |l3| by replacing c qi (t — r) and c q2 (t — r) 
with c qi (t) and c q2 (t) respectively, and setting the upper 
limit of integration over r to infinity ,14j. 

It is clear by inspection of Eqs. J3J and JSJ that the sys- 
tem dynamics can be interpreted in terms of four-wave 
mixing processes involving two matter-wave fields and 
two optical fields. For example, the dynamics of the 
atomic quasi-mode q results from the combined effects 
of the laser and a vacuum field mode together with an 
additional atomic quasi-mode q', while the dynamics of 
a scattered photon mode {k, A} results in turn from the 
combined effects of the laser field and two atomic modes. 

For optical pulses of short duration, we can assume 
that the condensate population is undepleted and replace 
the operators Co and cj by the c-number \/No, where 
Nq is the number of condensate atoms. Inserting then 
Eq. J5J into Eq. (QJ results in, 

c q - Aj(i)c q + S q (<)c f _ q + r q (t), (6) 

where q ^ and we have only kept the linear terms in- 
volving c q and c q , consistently with the undepleted pump 
approximation. Here 

Mt) = ^^El^,A| 2 [l n q( k )I^K-^+^) 
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and Bq(t) = — A-q exp(2iu q t). 

The Langevin noise operators r q (t), which account for 
the vacuum electromagnetic fluctuations through the ini- 
tial photon operators -BkA(O) and S kA (0) and are respon- 
sible for the initiation of the scattering process 2] , have 
the form r q (i) = exp(iw g i) [/ q (i) - /_ q (t)], where 

k(t) -*^^E-9kASkA(0) e ^— »%(k). 
k,A 

The coefficients Aq and Bq each contains two terms, 
corresponding to the four processes illustrated in Fig. ^ 
The two processes described by Aq are resonant, while 
those described by S q are characterized by an energy 
mismatch 2hio qi as evident from the expression for Bq. 
Fig- n shows that Bq describes a pair production process 
where two condensate atoms are scattered into the quasi- 
modes q and — q. Such a process is known in quantum 
optics to result in the production of entangled or squeezed 
particle pairs. This process, which leads to the genera- 
tion of backward-scattered atoms, will clearly occur for 
interaction times short enough that energy-momentum 
conservation (or phase-matching) is not yet established. 
This explains why back-scattering was not observed in 
the quasi-cw experiments of Ref. .lj. Furthermore, Aq(t) 
is positive (negative) for forward- (backward-)scattered 
atoms. Hence, the onset of backward scattering relies 
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fully on the pair production characterized by Bq(t). 
Equation JSJ readily gives the quasi-mode populations 



(c q c,) as 
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where the appearance of the anomalous density rn q = 
77i_ q = (c_ q c q ) is a clear signature of the quantum cor- 
relations between the modes q and — q, and the quan- 
tities Afq and M.q result from the noise operators r q . 
Their explicit forms can be obtained using the quantum 
regression theorem as 



FIG. 1: Plots Ai and A2 {B\ and B2) are the two scattering 
processes represented by the coefficient Aq (S q ). "0" and 
"±q" label the atomic momentum states with "0" being the 
condensate mode. The solid and dashed arrows represent the 
pump and scattered photon, respectively. 
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and Mq = -(A/" q +AL q ) exp(i2ui q t)/2. 

For completeness, we also give the expression for the 
number of scattered photons in mode {k, A}, 
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Equations and (JSJ form a closed set that can be 
solved numerically. The atomic recoil frequencies to q 
are many orders of magnitude smaller than the laser 
frequency u>l, so that conservation of energy requires 
that the wave numbers of the scattered photons satisfy 
k ~ fci. Furthermore, for the large condensates consid- 
ered here, the slowly-varying amplitude approximation 
mentioned earlier implies that n q (k) is only significantly 
different from zero for k + q kj, ~ 0, indicating momen- 
tum conservation for the scattering process. Under such 
conditions, the atoms are scattered from the condensate 
in two dipole emission halos, represented schematically 
by the two circles of radius k L in Fig. [21 

We evaluate the coefficients in Eqs 
by following the procedure of Ref. 
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12 fcfJd| 2 /(47r?iA 2 ) characterizes the strength of the 
atom-photon interaction, and 
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Here, w and Z are the dimensions of the cylindrically sym- 
metric condensate along the radial and axial directions, 
respectively, and 9^ is the angle between k(= k^ — q) 
and the long axis of the condensate. For a cigar-shaped 
condensate, I 3> w, the scattered photons are predom- 



scattered atoms 




FIG. 2: Schematic of the dipole emission pattern. The circle 
with solid (dashed) line is traced out by the momenta q of 
the forward- (backward-) scattered atoms. 



inantly along that axis, Dicke's end-fire mode, so that 
most of the scattered atoms gain momenta at roughly 
45° from the long axis of the condensate 0, 0, |n| . 

Figure shows the ratio of the backward- and 
forward-scattered atom numbers along this 45° angle. 
This ratio is always less than 1 since the backward scat- 
tering process has an energy mismatch AE = 2hu q while 
the forward process is resonant 11]. Fig. 03 shows 
the angular distribution of the scattered atoms. The 
forward-scattered atoms form an almost perfect symmet- 
ric distribution about the 45° angle. This must be con- 
trasted to the backward-scattered atoms, whose distri- 
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FIG. 3: a) Ratio of the backward- and forward-scattered atom 
numbers at angle <j> = 45° as a function of time, in units of 
1/iOr, where uj r = hk 2 L /{2M). The lines from top to bottom 
correspond to 7 = 500, 250, 150 and 50cj r . b) Angular dis- 
tribution in arbitrary units of the forward- (solid line) and 
backward-scattered (dashed line) atoms. Here 7 — 50u> r and 
t — 0.7 /lo t . The results are obtained by numerically integrat- 
ing Eqs. and JSJ with initial conditions n q = m q = using 
a 4 th -order Runge-Kutta method. 

bution shows an asymmetry, with more atoms scattered 
into 4> > 45° than into <j> < 45° (see Fig. Backward 
scattering favors larger angles because the energy mis- 
match is smaller for 4> > 45° than that for <f> < 45°, 
since AE = 2hui q oc cos 2 <f>. A similar asymmetry in the 
behavior of the forward- and backward-scattered atoms 
was observed experimentally |Tlj . and was attributed to 
the fact that the two scattering processes occur mainly 
at different locations along the condensate. Because the 
Markov approximation neglects the effects of retardation 
in the description of the end-fire mode dynamics [l6j| . 
that mechanism is absent from our analysis. In practice, 



both phase matching and spatial effects are probably at 
play, and they could in principle be distinguished by vary- 
ing the time of flight before detection [l2j . 

In summary, we have presented a theoretical investiga- 
tion of a zero-temperature condensate interacting with a 
far off-resonant laser pulse. Our work is valid for inter- 
action times short enough that the condensate remains 
undepleted and one can retain only linear terms involv- 
ing atomic modes of finite momenta. Our analysis shows 
that while the resonant Raman scattering illustrated in 
the plots Ai t 2 of Fig.^is the dominant generation mecha- 
nism of forward-scattered atoms, the backward-scattered 
atoms result from the quantum-correlated pair produc- 
tion processes of plots B\ i in Fig. 1S| . The lack of 
phase-matching for such pair production processes re- 
sults in distinct angular distributions for the backward- 
and forward-scattered atomic peaks, in good qualitative 
agreement with the experiment of Ref. . 

Schncble et at interpret the existence of backward- 
scattered atoms as resulting from Kapitza-Dirac diffrac- 
tion of matter waves off the optical gratings formed by 
the pump laser and the end-fire modes [llj. While this 
point-of-view is consistent with our correlated pair pro- 
duction picture, we emphasize that the strength of the 
optical grating depends on the intensity of the end-fire 
modes, which in turn depends on the population of the 
atomic quasi-modes [see Eq. (JjJJ]. It is this interplay of 
the optical and atomic fields that renders the backward- 
and forward-scattered atoms correlated. It also shows 
that the two diffraction processes of the system — optical 
diffraction off the atomic grating and atomic diffrac- 
tion off the optical grating — are coherently mixed and 
hence inseparable. For the same reason, the two counter- 
propagating optical end-fire modes must also exhibit 
quantum correlations. 
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